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1.  DISCUSSION 

V 

Two  approaches  to  finite  element  analysis  are  now  widely  recognized  in  the 

O'"'  u 

engineering  and  mathematical  communities.  In  both  approaches  the  domain  K  is 
divided  into  simple  convex  subdomains  (usually  triangles  or  rectangles  In  two 
dimensions,  and  tetrahedra  or  bricks  in  three  dimensions)  and  over  each  sub¬ 
domain  the  unknown  (displacement  field)  is  approximated  by  a  (local)  basis 
function  (usually  a  polynomial  of  degree  p) .  Basis  functions  are  required  to 
join  continuously  at  boundaries  of  the  subdomains  In  the  case  of  planar  or  3 
dimensional  elasticity,  or  smoothly  in  the  case  of  plate  bending.  The  differ¬ 


ence  between  the  two  approaches  lies  in  the  manner  in  which  convergence  is 
achieved.  ^These  two  approaches  are:  )  — 

1.  The  h-version  of  the  finite  element  method.  In  this  approach  the  degree  p 
of  the  approximating  polynomial  is  kept  fixed,  usually  at  some  low  number 
such  as  2  or  3.  Convergence  is  achieved  by  allowing  h,  the  maximum  diameter 
of  the  convex  subdomains,  to  go  to  zero.  Estimates  for  the  error  in  energy 
have  long  been  known.  In  all  of  these  estimates  p  Is  assumed  to  be  fixed  and 
the  error  estimate  Is  asymptotic  in  h,  as  b  goes  to  zero. 

2.  The  p-verslon  of  the  finite  element  method.  In  this  approach  the  subdivi¬ 
sion  of  the  domain  fl  is  kept  fixed  but  p  is  allowed  to  Increase  until  a 
desired  accuracy  is  attained.  The  p-version  is  reminiscent  of  the  Rltz 
method  for  solving  partial  differential  equations  but  with  a  crucial  distinc¬ 
tion  between  the  two  methods.  In  the  Ritz  method  a  single  polynomial 
approximation  is  used  over  the  entire  domain  ft  (ft,  in  general,  is  not  con¬ 
vex)  .  In  the  p-version  of  the  finite  element  method  polynomials  are  used 
as  approximations  over  convex  subdomains.  This  critical  difference  gives 
the  p-version  a  much  more  rapid  rate  of  convergence  than  either  the  Ritz 


method  or  the  h-version . 


The  p-verslon  of  the  finite  element  method  requires  families  of  polyno¬ 
mials  of  arbitrary  degree  p  defined  over  different  geometric  shapes.  Polyno¬ 
mials  defined  over  neighboring  elements  join  either  continuously  (are  in  C°) 
for  planar  or  three  dimensional  elasticity,  and  smoothly  (are  in  C^)  for  plate 
bending.  In  order  to  Implement  the  p-version  efficiently  on  the  computer, 
these  families  should  have  the. property  that  computations  performed  for  an 
approximation  of  degree  p  are  re-usable  for  computations  performed  for  the 
next  approximation  of  degree  p  +  1.  We  call  families  possessing  this  property 
hierarchic  families  of  finite  elements. 

The  h-version  of  the  finite  element  method  has  been  the  subject  of  inten¬ 
sive  study  since  the  early  1950's  and  perhaps  even  earlier.  Study  of  the 
p-version  of  the  finite  element  method,  on  the  other  hand,  began  at  Washington 
University  in  St.  Louis  in  the  early  1970's.  Research  in  the  p-version  (formerly 
called  The  Constraint  Method)  has  been  supported  in  part  of  the  Air  Force  Office 

of  Scientific  Research. 

Recent  Research  Accomplishments  are: 

1.  The  proof  a  theorem  which  establishes  that  for  plate  bending  problems  (in  which 

-continuity  is  required)  the  rate  of  convergence  of  the  p-version  is  twice 
that  of  the  h-version.  This  supplements  an  earlier  theorem  which  was  proved 
for  the  C°  case. 

2.  Hierarchic  families  of  C  -triangular  elements  have  been  implemented  on  the 
computer  using  pre-computed  arrays. 

3.  The  difficult  problem  of  the  finite  element  analysis  of  a  rhombic  plate  has 
been  solved  using  both  the  p-version  and  the  h-version. 

Work  has  begun  on  an  efficient  post-processing  technique  to  compute  accurate 
displacements  and  stresses,  and  to  calculate  stress  intensity  factors. 
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